Quasi-particle Tunneling Through a Barrier in the Fractional Quantum Hall Regime 
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Tunneling of fractionally charged quasi-particles (QPs) through a barrier is considered in the 
context of a multiply connected geometry. In this geometry global constraints do not prohibit such 
a tunneling process. The tunneling amplitude is evaluated and the crossover from mesoscopic QP- 
dominated to electron-dominated tunneling as the system's size is increased is found. The presence 
of disorder enhances both electron and QP tunneling rates. 
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One of the most remarkable facts about the fractional 
quantum Hall effect (FQHE) is the existence of fraction- 
ally charged quasi-particles (QPs). Their dynamics is 
manifest in a host of physical phenomena, whose ob- 
servation strongly supported the veracity of Laughlin's 
theory [1]. It has been pointed out [2, 3] that QP tunnel- 
ing is distinctly different from electron tunneling. Pertur- 
bative renormalization-group analysis [4] has indicated 
that in the weak backscattering limit inter-edge tunneling 
through the FQHE liquid is dominated by QP tunneling. 
These predictions have been confirmed by experiments 
[5] . In the opposite limit of strong backscattering (nearly 
disconnected FQHE systems coupled by weak tunneling 
through an insulator), the same RG analysis would have 
predicted that tunneling should be dominated again by 
QP tunneling. Common wisdom, however, has it that in 
this limit only electron tunneling is possible. The ratio- 
nale for that goes as follows: consider two FQHE puddles 
weakly connected through tunneling. The total number 
of electrons on each puddle is (nearly) a good quantum 
number; hence it must be an integer. QP tunneling would 
render this number non-integer, therefore such a process 
must be excluded. 

Our starting point here is to note that there are se- 
tups where the above mentioned "global constraint" (i.e. 
the number of electrons on each side of the barrier be- 
ing an integer) does not exclude a-priori QP tunneling 
through a potential barrier. The common wisdom al- 
luded to above needs then to be re-examined. Studying 
these setups is particulary interesting in view of recent 
experimental results [6] which suggest the coexistence of 
both electron and QP tunneling under strong backscat- 
tering conditions. 

Consider first the annulus depicted in Fig. 1(a). 
Clearly, the passage of a QP through the barrier would 
not violate the 'global constraint'. There are two pos- 
sible trajectories (for the presumed noiseless incoming 
current) to traverse the system: either following the edge 
adiabatically, or by tunneling through the barrier. The 
outgoing current would then be noisy [7]. By analyzing 
this non-equilibrium noise one may detect the effective 



charge involved. 

Our extensive analysis, performed on a torus geometry 
described below, leads to three main results: (i) For our 
multiply-connected geometry and in the presence of a real 
potential barrier we confirm the existence of QP tunnel- 
ing, which decreases rapidly with system size (Eq. 3). (ii) 
We study the amplitude of such QP tunneling processes 
and identify the crossover, in terms of the system's pa- 
rameters, between the electron-tunneling-dominated and 
the QP-dominated regimes, (iii) We show that in FQHE 
systems, the presence of impurities may enhance both 
electron and QP tunneling amplitudes (Eq. 5), in the 
spirit of the Shklovskiii-Li-Thouless mechanism [8]. 
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FIG. 1: (a) Proposed experimental setup for an annulus, allow- 
ing the tunneling of both QPs and electrons through a barrier. 
The edge states are marked, (b) The torus geometry studied in 
this work; x and y represent two Cartesian coordinates, the unit 
cell being 2nR X L\ $i and $2 are two Aharonov-Bohm fluxes 
threading the torus. A Hall liquid (of density 1/m) covers the 
torus everywhere, except around the barrier which is initially dry 
(" extended- holes" ) . Here the barrier potential is circularly sym- 
metric. 

The torus geometry. To facilitate our analytical study 
and to avoid complications emanating from the sys- 
tem being open, we hereafter focus on a setup defined 
on a torus, spanned by the two periodic coordinates 
< x < 2ttR and < y < L (Fig. 1(b)), with a uni- 
form magnetic field perpendicular to the surface of the 
torus. On top of the torus surface we introduce a cir- 
cular potential ridge (Vo). The tunneling investigated 
here is between the two sides of this barrier. The main 
steps in the analysis are: (i) We first construct modified 
Laughlin-Haldane-Rezayi states which corresponds to a 
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bulk filling factor 1/m (the "wet area") and a "dry area" 
(made of extended holes) where the electron density is 
suppressed, ideally to zero. The ground state configura- 
tion is obtained by maximizing the overlap of the dry area 
of the many-body configuration with the barrier [9]. (ii) 
The Hamiltonian, hence the wave- functions, depend on 
two gauge fluxes, $1 and $2 (Fig. 1(b)). By adiabatically 
increasing $1 any many-body configuration will slide 
rigidly in the y-direction, giving rise to a change in its 
energy. As the levels of two different configurations cross, 
the ground state of the system changes abruptly. Below 
we show that such a change corresponds to a tunneling 
event. The set of many-body energy levels is plotted in 
Fig. 2(a). (iii) To enable tunneling we break the circular 
symmetry (in the x-direction) , introducing an additional 
asymmetric potential (Vi). This gives rise to finite ma- 
trix elements between different configurations. Avoided- 
crossing gaps in the energy-flux spectrum (Fig. 2(a)) are 
a manifestation of tunneling: the period in flux reflects 
the nature of the particle that tunnels. Below we calcu- 
late these tunneling matrix elements, (iv) We demon- 
strate quantitatively how the presence of a multitude of 
(5-function impurities enhances the tunneling. 




FIG. 2: (a) Intersecting energy curves as a lunction ol the flux <t>i. 
Avoided crossings correspond to tunneling of electrons and QPs as 
indicated. For QP (electron) tunneling the periodicity of the adia- 
batically varied ground-state energy is (f>o (3</>o) [3]. (b) Schematic 
density profile of the initial ground state (solid line) and the 
first excited state \Eq (dashed line). As <E>i increases, these density 
profiles slide to the left, eventually becomes the ground state, 
(inset) An actual density profile for m = 3, N = 6, N h = 1. Densi- 
ties corresponding to other wavefunctions are rigid shifts thereof. 

Our initial Hamiltonian is H — Hq + Vint + Vq, where 
Hq includes the kinetic part as well as the magnetic field 
and fluxes and J7 int is the two-particle interaction. The 
barrier potential Vo is assumed to be sufficiently weak to 
exclude mixing with higher Landau levels. 

For Vq = and no dry area Haldane and Rezayi [10] 
have found a set of m degenerate Laughlin wave func- 
tions. The solution is obtained for a magnetic field which 
is quantized according to the Dirac's condition RL = 
£ H 2 Nrp, where N$ is the number of magnetic flux quanta 
perpendicular to the torus surface and £ H = yjTic/eB is 
the magnetic length. These wavefunctions are eigenfunc- 
tions of total quasi-momentum (TQM) [11] and, simi- 
larly to the Laughlin wave function [12], are exact (zero 



energy) ground states of the Hamiltonian with hard core 
interaction (i.e. V 2 <5(r)). Their ground state electron 
density is nearly uniform. Here — mN, where N is 
the number of electrons. 

Extended Hole Wavefunctions. To render the barrier 
(and its close vicinity) "dry", we tune the magnetic field 
to allow for Nh holes: A^ = mN + Nh- The lowest Lan- 
dau level consists of N^ single particle states \n > (cf. 
Ref. [10]). The density profile of each single particle state 
is approximately a Gaussian in the y-direction and uni- 
form in the x-direction. The distance between the guid- 
ing centers of adjacent states is L/N^. In the subspace of 
the lowest Landau level (spanned by ( N * ) possible Slater 
Determinants) the ground state is determined solely by 
the interaction. Diagonalizing the hard-core interaction 
results in a set of zero-energy ground states, each having 
Nh holes extended in the x-direction. As an example, 
consider Nh = 1- When the interaction term is diago- 
nalized one obtains N^ Laughlin-like states of zero en- 
ergy. Each of these states corresponds to a non-uniform 
electron density: the filling is 1/m almost everywhere, 
but the occupation of one of the single-particle quasi- 
momentum states \n > is suppressed: the area around 
this guiding center is "dry" (Fig. 2(b)). We denote such 
a many-body state with an extended hole at \n > by \t n . 
It is an eigenstate of the TQM. By sliding all guiding 
centers rigidly by 1 we increase the TQM by N, shifting 
— > ^ n +i. The above procedure is readily generalized 
to Nh > 1. Choosing the Nh holes to be contiguous leads 
to a dry area of a desired width [14]. 

tunneling: electrons versus QPs. The initial ground 
state whose dry area coincides with the barrier is denoted 
by 4"o- As the flux $1 is varied adiabatically this many- 
body configuration slides rigidly around the torus. The 
dry area of ^0 moves and its energy Eq($>i) increases. 
The dry area of whose TQM differs from that of ^0 
by N mod(A«ji), slides towards the barrier and its energy 
Ei($i) decreases. When $1 is increased by </>o/2, E ($i) 
and Ei($i) intersect, and the ground state of the system 
switches ^0 - * ^l- This corresponds to a shift of each 
single electron state byl: |n>^|n+l>. Since the 
average occupation of \n > is 1/m this process describes 
a shift of charge of e/m from one side of the barrier to the 
other, i.e. QP tunneling. The level-crossing degeneracy 
is lifted by breaking the circular symmetry of the poten- 
tial Vq: H — ► H + V\. The QP tunneling matrix element 
is < ^0 1 Vi I \I/ 1 >. By analogy, the two many-body states 
*o and \l/ m differing in their TQM by iJV mod(A^) will 
cross when $1 increases by m^o/2, hence < ^ol^il^m > 
is the matrix element for electron tunneling. 
In order to evaluate these tunneling matrix elements 
we use the fact that V\ is a single-particle potential. 
The overlap < ^o|Vi > (p = l,m for QP, elec- 
tron tunneling respectively) consists of contributions 
from the respective Slater determinants components 
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\ki, . . ■ , fejv >G *o and \£i, . . . . 
tical except for a single pair k. 



In >£ * P which are iden- 
I. The difference in TQM 



l-k = pJVmod^). 



(1) 



Taking VI = Vi<5(x) renders the procedure particularly 
simple [15]. Then < fc|Vi|^ >= v p depends only on the 
difference £— k. The tunneling amplitude can then be cast 



as % =< #n|Vi 



* P >= 



fp(L)v P 



where f p is a combina- 



torial factor resulting from summation over all possible 
pairs satisfying (1), and is calculated numerically. For v p 
one readily obtains 
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(2) 



To leading order v p oc exp [-{ P N mod/V This 

reflects the overlap of two Gaussians separated by a dis- 
tance pNL/N^ on the torus (Fig. 3(a)). For a QP the sep- 
aration is of the order L/m (N^ = mN + Nh), so v p scales 
as exp[— (L/2m£ H ) 2 ]. By contrast, for electrons the sep- 
aration distance (defined modulo L) is N^L/N^,, hence 
v p scales as exp[-{LN h /2N^t a ) 2 ] = exp[-(N h £ H /R) 2 ], 
which is L- independent (where we used the Dirac con- 
dition). Likewise we find that the factor f p is roughly 
system-size independent for electron tunneling, while it 
rapidly decreases (Gaussian-like) for QPs. 
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FIG. 3: Tunneling involves initial and final states whose TQM 
differ by pN, p = 1,3 for QPs and electrons respectively. Here 
m = 3. The respective matrix elements involve the overlap of two 
single-particle states (Gaussians), whose separation is pNL/N^. 
For QPs the latter is of the order of L/m (solid arrow, black and 
white Gaussians). Thus the overlap factor for QPs, v p =i, strongly 
decreases with L. By contrast, for electrons (dashed arrow, black 
and grey) the distance (defined modulo L) is N^L/N^ << L and 
the dependence on L is negligible, (b) The calculated tunneling 
probabilities for electrons (73) and QPs (Ti), for m = 3, = 1. 
T3 is practically independent of system size (squares), while for QPs 



(circles). The solid line is a fit with a pa 0.07. 



We thus summarize 



QPs (p = 1) 



(3) 



L independant, electrons (p = 3) 
For Nh — 1 (and N = 2, . . . , 6) we obtain numerically 



(= l/4m) in the thermodynamic limit (cf. Eq. (5) for 
Ni m p = 1). Note that the factor f p= i further suppresses 
the single-particle term e - {L/2al H )2 -> e -( L / 2e H ) 2 / m . 

impurity-assisted tunneling. We next introduce im- 
purities into the system, Ti. = Hq + Umt + V%, where 



2^.j=i 



ZQj) represents Ni mp localized 



impurities (Vi mp > 0). Shklovskii and Li & Thouless 
[8] have shown that for non-interacting electrons, impuri- 
ties modify the Gaussian decay of the edge-to-edge Green 
function into exponential. This result is generalized here 
to the FQH regime. To simplify the analysis we em- 
ploy the fact that QP tunneling may be interpreted as a 
non-local process taking place through the liquid, while 
electron tunneling takes place through both the poten- 
tial barrier and the liquid. Thus, the torus can be ef- 
fectively replaced by two cylinders whose circumference 
is 2ttR and whose lengths are Lf, arrier = N^L/N^ and 
Liiquid = L - L bamer = mNL/N^. For a cylinder with 
impurities the ground state wavefunction ^{ Zo } contains 
Nimp localized holes at the impurity positions 



N h N 

nn^ 

j=i k=i 



-iz k /R 



(4) 



I>?/ 2 i 



is Laughlin's 



(*£ = n^e-^-e-'^r e 
cylinder wave function). 

To obtain the various tunneling matrix elements we 
calculate the overlap between ^{zo} an d its shifted ver- 

sion * {zo} - n^ie-^/^zo}, V = l,m for QP and 
electron tunneling respectively. For Ni mp = 1 one recov- 
ers (up to prefactors) the perturbative result [16]: QPs 
tunnel more efficiently than electrons along a cylinder. 
We have evaluated this overlap numerically for systems 
with TV < 6 (for electrons tunneling at integer filling we 
have considered N < 17) and N imp < 4, for (i) impuri- 
ties equally spaced on a line and (ii) at random positions 
throughout the cylinder, averaging over ~ 1000 realiza- 
tions. For case (i) we find that 



<^{zo}l*{z } >' 



D -L 2 /12N it 
,-L 2 /iN im 



QPs 

electrons 



(5) 



agrees with the numerics. For case (ii) the decay fac- 
tor of the exponent is modified, but not the parametric 
dependence. We find that when the longitudinal impu- 
rity density A = L/Ni mp is held constant, the typical 
hopping distance is kept fixed, a Gaussian-to-exponential 
crossover takes place. This crossover can be understood 
in terms of multiple impurity-assisted tunneling. For a 
QP, as an example, e ~ L / 12N im. p ^ e ~ AL / 12 . 

QP-electron crossover. Studying this crossover is now 
experimentally feasible [17]. Here we present a frame- 
work to study it theoretically in a multiply connected 
geometry, e.g. the torus. As Lfearrter is varied (com- 
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a w 0.07 (Fig. 3(b)). We expect this to approach 1/12 pared with Lu qu id) the QP tunneling 7{ competes with 
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the electron tunneling By calculating T\ and T p , as 
explained above, varying the number of particles (modifi- 
cation of Lii qui( i) and the number of extended holes (mod- 
ification of Lt, arrier ), we obtain the ratio of the electron- 
tunneling amplitude to that of QP tunneling (Fig. 4). 
This allows us to determine the separation between 
the electron-tunneling dominated and the QP-tunncling 
dominated regimes (dotted line). This can be compared 
to the solid curve obtained by taking the estimates T\ ~ 
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FIG. 4: Electron-QP crossover: Plotted is the interpolation 
of In (T3/T1) for discrete values of N and Nh corresponding to 
Lu q uid/ m an d Lf >arrier respectively (measured in units of £ H 2 /R\ 
here R = \/l~47r). Dashed line: the crossover curve T3 = T\. 
Solid line: an approximated crossover-curve obtained by taking 
-£?. -J121 2 -l 2 IU 2 , -L 2 . JU 2 

e liquid 1 H — g barrier' H 4- e liquid' H 

discussion. Keeping the barrier size fixed and increas- 
ing the torus length we have found that the QP tunnel- 
ing amplitude decreases while the electron amplitude is 
mostly unaffected. This supports the picture that the 
QP tunneling through a barrier [18], while in principle 
possible, is a mesoscopic effect. It may be interpreted 
as a QP leaping through the liquid around the barrier. 
In the thermodynamic limit the QP tunneling amplitude 
vanishes, in accordance with common wisdom. The de- 
pendence of the tunneling amplitude on length scales 
does not conform to the scaling resulting from the RG 
treatment [4]. It is strongly modified by the multiple 
connectedness of the system. Adding disorder enhances 
the tunneling amplitudes. As can be seen from Eq. (5), 
special arrangements of the impurities can lead to even 
stronger enhancement (e.g. increase the linear density of 
Ni m p, while keeping the two-dimensional density fixed). 
We believe that the best candidate to test the ideas out- 
lined here is the annulus geometry (Fig. 1(a)). The rele- 
vance of the current results to the annulus geometry will 
be explored in future studies. 
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